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Abstract 

We investigate the possibility of calculating the parameters of few-body resonances using 
the oscillator trap boundary conditions. We place the few-body system in an oscillator trap 
and calculate the energy spectrum and the strength function of a suitably chosen transition. 
Broader resonances are identified as Lorentzian peaks in the strength function. Narrower 
resonances are identified through the pattern of avoided crossings in the spectrum of the 
system as function of the trap size. As an example we calculate and 0;^ resonances in 
"^^C within the 3a model. 

1 Introduction 

Few-body resonances are often calculated using the complex scaling method where the reso- 
nances are identified as generalized complex eigenvalues of the Hamiltonian with the corre- 
sponding generalized eigenfunctions (see e.g. |il]l2| and references therein). The method has the 
advantage of having a simple boundary condition: the few-body wave-function vanishes at large 
distances. However, it also has certain disadvantages. Complex arithmetics and algorithms are 
generally slower and complex matrices need more computer memory. Calculating extremely 
narrow resonances is difficult as it demands calculations of the eigenvalues with exceedingly 
high accuracy. Interpretation of the generalized eigenfunctions is also not trivial P^, especially 
for heavy complex scaling needed for broader resonances. 

In this contribution we investigate the possibility of calculating the parameters of few-body 
resonances using the same simple boundary condition as in complex scaling method but working 
with real energies and real wave- functions. 

We place the few-body system in an artificial oscillator trap of length h which is significantly 
larger than the characteristic length of the few-body system. We then calculate the (discrete) 
spectrum of the system in the trap and estimate the strength function of a certain transition 
from a suitable chosen initial state to the positive-energy states of the system in the trap. 

The broader resonances with width F > /i^/(2m5^), where m is the characteristic mass of 
the few-body system, can be identified as Lorentzian peaks in the strength function. 

The narrow resonances with F < h?/{2mb'^) need an investigation of the spectrum of the 
system in the trap as function of the trap length. When an energy level in the trap, following 
its general behavior as approaches the system's resonance level to within its width, the 
two levels interfere and avoid crossing. The pattern of avoided crossings in the spectrum of the 
system in the trap reveals the position of the narrow resonance. The width of the resonance 
can be estimated from the size of the region of avoided crossing, or, more precisely, from the 
variation of the energy levels with respect to the trap length. 
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The approach is similar to the box method (also called the stabilization method) [31 d]. 
However, the difference is that we use an oscillator trap instead of a box and that we resort 
to strength function method for broader resonances where the avoided crossings method is less 
reliable. The box boundary condition is more complicated as the wave-function has to vanish 
identically at the box boundary which for few-body systems is a multi-dimensional surface. 
The oscillator trap can be potentially used in stochastic variational calculations with correlated 
Gaussians [5j. 

As an example we apply the approach to the 3a system in the J'^=0^ channel where there 
exist a narrow, 0^, and a broader, Oj^, resonance. We show that the approach allows to reliably 
calculate the two resonances in this system. 



2 The few-body system and the trap 

We consider the 3a system with the total angular momentum and parity J'"=0^. The Ali- 
Bodmer type a-a potential is taken from [6], 



Vc^air) = (l25P,=o + 20Pi=2)e-(T^)' -30.18 e 
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where all energies are in MeV, all lengths in fm. Pi is the projection operator onto a state with 
relative orbital momentum /, and r is the distance between a-particles. In addition a three-body 
force 

^3(^5) = -76MeVexp(-pV(4fm)2) , (2) 

is employed to simulate the contribution of "compound nucleus" degrees of freedom at shorter 
distances where all three a-particles overlap. The three-body force is defined in terms of the 
hyper-radius 

p^ =— (3) 

where rj are the cm. coordinates of the a-particles, m=939 MeV is the chosen mass scale and 
ma=3.97m. 

The system is placed in an oscillator trap 

where the trap length h is varied around 30-40 fm. 

The three-body problem in the trap is solved using the adiabatic hyper-spherical method 
(see e.g. [2J and references therein). First for every fixed hyper-radius p the eigenvalue problem 
for all remaining variables (denoted collectively as hyper-angles VC) is solved and the spectrum of 
hyper-angular eigenvalues ei(/?) together with the angular eigenfunctions <^j(/9, f]) are obtained. 
The functions '&i(/o, f^) are then used as a full basis in the space and the total wave- function 
^ is represented as a series 

00 

V^(p,f]) = ^/;,(p)$,(p,Q) , (5) 

j=i 

where the expansion coefficients /i(p) are obtained by solving the hyper-radial equations where 
the eigenvalues ei(p) serve as effective potentials. 
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3 Strength function 



A resonance can be identified as a peak in a reaction cross-section with approximately Lorentzian 
shape. The amphtude of a quantum transition, caused by an operator F, from some initial state 
ipa into one of the discrete state ipn of the system in the trap, is given in the Born approximation 

as 

Mn^a = {4^n\F\tPa) ■ (6) 

Since a resonance per definition must be seen in any reaction channel, the particular choice of 
the excitation operator F and the initial state ipa should be irrelevant as soon as the matrix 
element does not vanish identically. We thus choose the initial state in the form of the large p 
asymptotics of a bound three-body state [Till], 

f,{p) = p-'/'eM-p/b3), (7) 

in every hyper-radial channel i. The constant 63=4 fm is chosen close to the size of the bound 
state of 3a system. The excitation operator is taken as 

F = p\ (8) 
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Figure 1: Left: the strength function S{E) as function of the 3a energy E for different oscillator 
lengths b. Right: the peak at 4.3 MeV is fitted with a Lorentzian, Const/[(E' - Er)"^ + T^/^, 
where Er=4:.3 MeV and r=0.9 MeV. 

The cross-section of a reaction into the final states with energies E it is determined by 
the strength function, defined as 

^(^) = ^ E \Mn^a\'. (9) 

The energy bin size AE has to be chosen on the one hand small enough as not to smear out 
the essential features of the cross-section, and on the other hand large enough to include many 
states. In our calculations the energy bins include four states each. 

The calculated strength function is shown on Figure [TJ It reveals a broader peak at 4.3 MeV 
and a narrow unresolved peak at ~0.4 MeV. In the region of the broader peak the strength 
function is well converged with respect to the trap length, and the bin size is quite appropriate 
for the description of the width of the peak as there are many points within the peak region. 

The shape of the peak is well described by a Lorentzian 

S{E) ^ (10) 

{E-Err + ^ 



3 




b, fm b, fm 



Figure 2: Left: the spectrum of the 3a system in an oscillator trap as function of the trap size 
b in the region of the narrow peak on Figure (T) the sequence of avoided crossings indicates a 
resonance at ~ 0.38 MeV; Right: zoom-in into the region of avoided crossings: the resonance 
energy is fitted with E{b) = Er + Kb'"^ with Er = 0.38435 MeV and i^=7.471^ MeV fm^. 



with i5^r=4.3 MeV and r=0.9 MeV. These numbers are consistent with [2]. 

The narrow peak at in the strength function ~0.4MeV is represented by only one point. The 
position of the point reveals the resonance energy but not the width. To resolve the width at 
least several points are needed within the peak region. For resonances width exceedingly narrow 
width r this would demand unreasonably large trap lengths of the order b ~ y^^y(^^j^ry_ 

However instead of the strength function the avoided crossings method can be used to 
calculate narrow resonances using reasonably sized traps. 

4 Avoided crossings 

For large trap lengths the energy levels in the trap scale with the trap size as ^ oc Varying 
the trap size a level in the trap can be moved close to the resonance level of the system. If 
the resonance were behind a completely impenetrable barrier (thus having a vanishing width) 
there would be no interference through the barrier between the resonance and the state in the 
external trap. The resonance would then be insensitive to the trap size. The spectrum of the 
system in the trap, as function of the trap size, would thus show the trap levels scaling as 
and crossing the resonance energy represented by a horizontal line. 

If the barrier has small but finite penetrability the trap level approaching the resonance to 
within its width becomes perturbed by the resonance resulting in the "repulsion" of the two 
interfering levels. This shows up as a sequence of avoided crossings in the spectrum of the 
system in the trap as function of the trap size. 

Indeed the spectrum of the 3a system in the trap reveals such a sequence of avoided crossings 
in the vicinity of the narrow resonance, see Figure [2] (left). 

The resonance state gets a contribution from the oscillator potential ^ which at large b 
is proportional to This contribution can be determined by a fit E'r + Kb~^ through the 
resonance energies as shown on Figure [2] (right). The fit also provides the asymptotic estimate 
of the resonance energy Er = 0.38435 MeV. 

The figure shows the reduced energy E'^ = E^ — Kb~^ where the oscillator contribution 
is subtracted. It is possible to estimate the width F of the resonance from the plot assuming 
that the energy region where the avoided crossing takes place is determined by the width of the 
resonance, 
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Figure 3: The reduced energy E'^ = E3 — Kb~^ of tlie tliird level of the 3a system in an oscillator 
trap as function of the trap length b in the region of the resonance Er = 0.38435 MeV (indicated 
as a horizontal line) where the parameters K and E^ are from the fit of the resonance energy 

arctan S^'t where b^, A^s, and T 
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on Figure [5J The energy E'^ is fitted with the curve 
are fitting parameters. The fit gives F = 77 eV. 



where Ab is the distance between the neighboring avoided crossings. 

However, instead of numerical differentiation it is better to estimate the width by fitting the 
calculated energies with the curve 

— arctan — — ^ , (12) 



Abn E'^ - Er 



where bn, Abn, and F are fitting parameters. Figure O (right) shows such a fit for E^. The fit 
gives F = 77 eV. This value is consistent with the estimates of 10-30 eV in ^ |9j taking into 
account that our three-body potential provides a slightly higher Er. 



5 Conclusion 

Using the two lowest J^=0^ resonances in the 3a system as an example, we have investigated 
the possibility of calculating the energies and widths of few-body resonances by placing the few- 
body system in an artificial oscillator trap. The oscillator trap has particularly simple boundary 
condition and can be potentially used in stochastic variational calculations with correlated 
Gaussians. 

We have shown that broader resonances with the width F > h? /{2mb'^), where b is the 
trap size, can be identified as Lorentzian peaks in the strength function of a suitably chosen 
"gedanken" transition. Narrower resonances can be identified through the pattern of avoided 
crossings in the spectrum of the system in the trap as function of the trap size. 
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